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Abstract 

We study universality properties of the Weingarten hyper-cubic random sur¬ 
faces. Since a long time the model of hypercubic random surfaces with a local 
restriction forbidding surface self-bendings was thought to be in a different uni¬ 
versality class from the unrestricted model defined on the full set of surfaces. 
In this paper we show that both models in fact belong to the same universality 
class with the entropy exponent 7 = 1/2 and differ by the finite size effects 
which are much more pronounced in the restricted model. 


The model of hyper-cubic surfaces was proposed by Weingarten Q as a 
non-perturbative regularization of the world sheet of bosonic strings. It was 
proven Q| that in general the string susceptibility exponent 7 is equal to 1/2 
for spherical surfaces, reflecting the fact that the ensemble is dominated 
by surfaces having the structure of branched polymers. It was also generally 
believed on grounds of the renormalization group arguments that one cannot 
change the exponent 7 by adding local terms to the action. Therefore it came 
as a surprise when results of the numerical simulations of a slightly modified 
model, where self-bendings surfaces were excluded from the partition function, 
gave 7=7 |@| in the case of a surface embedded in four dimensions. 

This result was subject of much speculation Q, because it provided the first 
non-trivial value, besides the mean field result 7 = i, in the series of positive 
7 discussed by Durhuus in ||^. On the other hand the result is not strictly in 
conflict with the analytic proof, as discussed in 

In this article we present new measurements of the exponent 7 . Since the 
earlier measurements signihcant progress has been made both in computer 
speed and in the development of numerical techniques for simulations of random 
geometries. Combining this allows us to go to systems two orders of magnitude 
larger than those used in the old measurements, and to gather very good statis¬ 
tics. Moreover we have at our disposal a very powerful method for extracting 
the exponent 7 based on the analysis of the baby-universe distribution. 

For this work we developed a new algorithm, which we will describe in a sep¬ 
arate article 0. This algorithm is capable to include self-bendings, although 
for most of the data presented in this letter they were suppressed by adding 
additional code. By reusing the 7 extraction method described in we show, 
that the results generated by our algorithm for surfaces of the size used in that 
work are in agreement with the results generated by their algorithm How¬ 
ever, as we will show our results on larger surfaces lead to completely different 
conclusions, namely that the exponent 7 is consistent with 1 / 2 , although with 
considerably stronger finite size corrections than in the case of unrestricted sur¬ 
faces. 

The model is defined by the partition function : 

= E = E (1) 

where the first sum runs over all hypercubic surfaces without self-bendings and 
fixed (spherical) topology, and the second runs over the number of plaquettes 
A, i.e. the surface area. The internal symmetry C(S) is important mainly for 
very small surfaces. 

In the large volume limit, the number of surfaces with a given area A is 
expected to behave like 

Ar(A) « (2) 

Our main interest is focused on the universal entropy exponent 7 , which also 
characterizes the fractal structure related to baby-universes. 
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To extract this exponent, the simulations are organized in a quasi-canonical 
way. The area of the system fluctuates but the measurements are taken at 
a fixed value Aq. We achieve this by adding an external quadratic potential 
— Aq)^ 0 to the exponent in the partition function]^. The optimal value of 
the parameter /3 is determined experimentally in an iterative procedure which 
after each run estimates an improved value of [3 from the A distribution. 

To measure 7 we use the baby universe method 0 . This method is based 
on the observation that the distribution of baby universes is independent of the 
critical value /Xc. By baby universe we mean a part of the surface which has the 
topology of a disc and is connected in such a way that the square of the disc 
perimeter I is much smaller than the disc area B : P B. When one looks 
at such a short loop from far enough away one sees it as a point to which two 
surfaces are pinned : one with area B and the other one, the complementary part 
of the whole surface, with area A — B. For large B and A — B the correlations 
of the parts of the surface on the two sides of the loop are negligible and one 
can approximate the baby universe distribution by the formula [0 : 

Na{B)^Mi{B)Ni{A-B) (3) 

where A/j (A) is number of discs with the loop I and area A. The number TVj (A) 
can be expressed by the number of spherical surfaces M{A). Take for example 
I as the perimeter of the square. We have then 

M{A) = {A + 1)M{A + 1). (4) 

This means that one can make such a disk by removing one square from the 
lattice. Generally, for short loops one then has 

M{A)^AM{A){l + o{l/A)}. (5) 

In our measurements we use loops of length two and four, being, respectively, 
the length of a double link and a square, i.e. the smallest possible objects. 
Using the approximate formula we get 

Na^B) - (BiA - ( 6 ) 

for the distribution of minbus with area i? in a canonical ensemble with a total 
area A. The exponential factor becomes a normalization factor independent 
of the minbu area B. Taking the logarithm of both sides and defining 

y = \n{B{A-B)fNA{B) 

X = \nB{A-B) ^ ’ 

one can extract 7 by htting to the linear equation y = jx + b. 

In fig(|^ we show a typical minbu distribution, where we plot y versus x. The 
total surface area is A = 10000. For small values of B the data points do not lie 
on a single line but wiggle around. This effect comes from the combination of 
the symmetry factor and the surface entropy. For large values of B one can see 
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Figure 1: The distribution of minbu-areas for A = 10000 


strong statistic noise which comes from the fact that the frequency of appearance 
of the large baby universes is very small. For example, for the surface A = 10000 
one has 7 10“"^ baby universes of area B = 4999 per surface. 

Due to finite size effects the curve is not quite a straight line. Therefore 
one actually measures an effective 7 e//j which should converge to 7 with larger 
areas A. When estimating the value of 7 one has to make decision which part 
of the distribution to use. We extract 7 from this distribution by looking at its 
behavior as a function of the lower cutoff, i.e. the smallest minbu-area used to 
fit 7 . As the value of 7 we take the plateau-value, as one can see for example 
in figure (^) for A = 10000. It is important to note that the plateau value sets 
in for baby universes larger than 1000 . 

We analyzed the minbu-distribution for different canonical volumes. In figure 
we show 7 e// as a function of the canonical volume A. For comparison 
we also performed simulations for the ensemble containing self-bendings. As 
discussed above, one can for this case show analytically that 7 = ^- For the 
model without self-bendings (lower curve) one clearly observes stronger finite 
size effects. The large area behavior of 7e//, however, is not affected. 

This observation is different from the results of the numerical experiment 
reported in |^. The reported value of 7 = 0.26 can be found in figure (||) for 
volumes A ss 200 which in fact correspond to the values used in that work. The 
method used there is based on the ratio of two shifted histograms of the area 
distribution f{A) for spherical surfaces obtained for different coupling constants 
in the action f3A. As before this eliminates the unknown parameter in the 
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Figure 2: The effective jeff as a function of the lower cut-off. 



Figure 3: The effective entropy exponent Jeff as a function of the canonical 
volume A for ensembles respectively including (upper curve) or excluding (lower 
curve) self-bendings. 
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exponent which disappears as it is independent of A. Thus one ends up with 
the formula : 

^ (l^) “ ^ 2 ))^ (8) 

where the histograms were measured for two different sizes A and A + 5. The 
method is is similar to the baby universe method in the sense that one gets rid of 
the inconvenient unknown exponential prefactor. Here, however, the price one 
has to pay is much higher because one has to use two independent histograms, 
which introduces a large statistical error. 

We have repeated measurements of 7 using this method and following pa¬ 
rameters : 

P A 

1.170 [170,280] 

1.160 [4890,5000] 

1.165 [4890,5000] 

1.170 [4890,5000] 

These are very similar to those used in the work [^. For each coupling we 
collect ten histograms, therefore we have 300 possibilities of calculating (||). The 
resulting distribution of 7 - values is is quite broad and turns out to react rather 
sensitive to a change of the surface-areas involved in the procedure. With the 
above parameters we get 'Jeff = 0.26(10), which is in very good agreement with 
the results presented in [^. This value is different from the branched polymer 
value 1/2 because 7 is measured at small volumes, where the finite size effects 
hide the real value, as can be clearly seen in figure (^). 

To summarize, we have measured the value of the string susceptibility ex¬ 
ponent 7 for hypercubic surfaces embedded in d = 4-dimensional target space 
for the ensemble of spherical surfaces with and without self-bending hinges. In 
both cases the results are compatible with the branched polymer value 7 = 1 / 2 . 
This means that the local infinite coupling to the external geometry, which 
suppresses self-bendings, is not sufficient to modify the fractal structure of the 
surface. These results do not rule out the possibility that one can influence 
the fractal structure with a generalized local coupling to the external geometry. 
Obviously with an infinite coupling to the external curvature the surface has to 
be (up to finite size effects) entirely flat, i.e. the Hausdorff-dimension is dh = 2 
in this case. This value is different from dh = 4, the value obtained for the 
model investigated in this work. Therefore the fractal structure is affected. It 
remains to find out, if a transition appears at finite couplings. 

We have developed a very general implementation of our new algorithm, 
which allows to study more general actions controlling self-bendings or other 
local terms in the action. Moreover, the program can be run already in the 
target space dimension two. We will present results of these simulations in a 
forthcoming publication. 
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